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Abstract 

We consider the class of non-Hamiltonian and dissipative statistical systems with dis- 
tributions that are determined by the Hamiltonian. The distributions are derived analyt- 
ically as stationary solutions of the Liouville equation for non-Hamiltonian systems. The 
class of non-Hamiltonian systems can be described by a non-holonomic (non-integrable) 
constraint: the velocity of the elementary phase volume change is directly proportional 
to the power of non-potential forces. The coefficient of this proportionality is determined 
by Hamiltonian. The constant temperature systems, canonical-dissipative systems, and 
Fermi-Bose classical systems are the special cases of this class of non-Hamiltonian systems. 
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1 Introduction 

The canonical distribution for the Hamiltonian systems was defined by Gibbs in the book 
"Elementary principles in statistical mechanics" [1], published in 1902. In general, classical 
systems are not Hamiltonian systems and the forces are the sum of potential and non-potential 
forces. Non-Hamiltonian and dissipative systems can have the same distributions as Hamil- 
tonian systems. The canonical distributions for the non-Hamiltonian and dissipative systems 
were considered in [2, 3, 4, 5, 7, 8, 6, 9, 10, 11]. 

The aim of this work is the extension of the statistical mechanics of conservative Hamiltonian 
systems to a wide class of non-Hamiltonian and dissipative systems. 

Let us point out non-Hamiltonian systems with distribution functions that are defined by 
^ ' the Hamiltonian. 

(1) In the papers [2, 3, 4, 5, 6, 9], the constant temperature systems with minimal Gaussian 
constraint are considered. These systems are the non-Hamiltonian systems that are described 
by the non-potential forces in the form f\ = — 7p« and the Gaussian non-holonomic constraint. 
Note that this constraint can be represented as an addition term to the non-potential force. 

(2) In the papers [12, 13], the canonical-dissipative systems are considered. These sys- 
tems are the non-Hamiltonian systems that are described by the non-potential forces F^ = 
—dG(H)/dpi, where G(H) is a function of Hamiltonian H. Note that the distribution functions 



1 



are derived as solutions the Fokker-Planck equation. It is known that Fokker-Planck equation 
can be derived from the Liouville equation [14]. 

(3) In the paper [15], the systems with non-holonomic constraint and non-potential forces 
= are considered. The equations of motion for this system are incorrect [16]. The correct 

form of the equations is derived in [15] by the limit r — > 0. This procedure removes the incorrect 
term of the equations. 

(4) In the paper [11], the canonical distribution is considered as a stationary solution of 
the Liouville equation for a wide class of non-Hamiltonian system. This class is defined by a 
very simple condition for the non-potential forces: the power of the non-potential forces must 
be directly proportional to the velocity of the Gibbs phase (elementary phase volume) change. 
This condition defines the general constant temperature systems. Note that the condition is a 
non-holonomic constraint. This constraint leads to the canonical distribution as a stationary 
solution of the Liouville equations. For the linear friction, we derived the constant temperature 
systems. The general form of the non-potential forces is not derived in [11]. 

(5) In the paper [17], the quantum non-Hamiltonian systems with pure stationary states 
are considered. The correspondent classical systems are not discussed. 

(6) In the paper [19], the non-Gaussian distributions are suggested for the non-Hamiltonian 
systems in the fractional phase space. Note that non-dissipative systems with the usual phase 
space are dissipative systems in the fractional phase space [19]. 

Khintchin [20] revealed the deep relation between the Gaussian central limit theorem and 
canonical Gibbs distribution. However, the Gaussian central limit theorem is non-unique. Levy 
and Khintchin have generalized the Gaussian central limit theorem to the case of summation 
of independent, identically distributed random variables which are described by long tailed 
distributions. In this case, non-Gaussian distributions replace the Gaussian in the generalized 
limit theorems. It is interesting to find statistical mechanics and thermodynamics that is based 
on non-Gaussian and non-canonical distributions [21, 22, 23, 24]. 

The aim of this paper is the description of non-Hamiltonian and dissipative systems with 
(canonical and non-canonical) distributions that are defined by Hamiltonian. This class can 
be described by the non-holonomic (non-integrable) constraint: the velocity of the elementary 
phase volume change must be directly proportional to the power of non-potential forces. The 
coefficient of this proportionality is determined by the Hamiltonian. These distributions can 
be derived analytically as solutions of the Liouville equation for non-Hamiltonian systems. The 
special constraint allows us to derive solutions for the system, even in far-from equilibrium 
states. This class of the non-Hamiltonian systems is characterized by the distribution functions 
that are determined by the Hamiltonian. The constant temperature systems [2, 3, 4, 5, 6, 9], 
the canonical-dissipative systems [12, 13], and the Fermi-Bose classical systems [13] are the 
special cases of suggested class of non-Hamiltonian systems. 

In section 2, the definitions of the non-Hamiltonian and dissipative systems, mathematical 
background and notations are considered. In section 3, we consider the condition for the non- 
potential forces. We formulate the proposition that allows us to answer the following question: 
Is this system a canonical non-Hamiltonian system? We derive the solution of the N-particle 
Liouville equation for the non-Hamiltonian systems with non-holonomic constraint. In section 
4, we consider the non-holonomic constraint for non-Hamiltonian systems. We formulate the 
proposition which allows us to derive the canonical non-Hamiltonian systems from the equations 
of non-Hamiltonian system motion. The non-Hamiltonian systems with the simple Hamiltonian 
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and the simple non-potential forces are considered. In section 5, we derive the class of non- 
Hamiltonian systems with canonical Gibbs distribution as a solution of the Liouville equation. 
In section 6, we consider the non-Gaussian distributions as solutions of the Liouville equations 
for the non-Hamiltonian systems. In section 7, we derive the analog of thermodynamics laws for 
the non-Hamiltonian systems with the distributions that are defined by Hamiltonian. Finally 
a short conclusion is given in section 8. 

2 Definitions of Non-Hamiltonian, Dissipative and Canon- 
ical Non-Hamiltonian Systems 

Let us consider the definitions of non-Hamiltonian and dissipative classical systems [25] , which 
are used for the formulation of our results. 

Usually a classical system is called a Hamiltonian system if the equations of motion are 
determined by Hamiltonian. The more consistent definition of the non-Hamiltonian system is 
connected with Helmholtz condition for the equation of motion. 

Definition 1. A classical system which is defined by the equations 

^ = G t , ^=F U (1) 
dt ' dt K ' 

where i = 1,...,N, is called Hamiltonian system if the right-hand sides of Eq. (1) satisfy the 
Helmholtz conditions for the phase space 

dd dG, dFi dG j dF t OF, 

dpj dpi dpj dqi c% c% 

Here Gi = Gi(q,p), Fj = Fi (q,p,a,t), where a is a set of external parameters. 

If the Helmholtz conditions are satisfied, then the equations of motion for the system (1) 
can be represented as canonical equations 

dqj _ dH_ djpj_ _ _dH_ . , 

dt dpi ' dt dqi ' 

which are completely characterized by the Hamiltonian H = H(q,p, a). In this case, the forces, 
which act on the particles are potential forces. 

If the functions Gi for the non-Hamiltonian system (1) are determined by the Hamiltonian 

Gi = j—, 4 
dpi 

and the Hamiltonian is a smooth function on the momentum space, then the first condition (2) 
is satisfied 

d 2 H d 2 H 
dpidpj dpjdpi 
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In this case, the second condition (2) has the form 

OF d 2 H 



dpj dq-idpj 



= 0. (5) 



In general, the second term does not vanish. For example, in the nonlinear one-dimensional 
sigma-model [27] the second term of the left-hand side of Eq. (5) is defined by the metric. 

Definition 2. A mechanical system is called non-Hamiltonian if at least one of conditions 
(2) is not satisfied. 

Let us consider the time evolution of the classical state which is defined by the distribution 
function pjy(g,p, a, t). The N-particle distribution function in the Hamilton picture (for the 
Euler variables) is normalized by the condition 

J PN(q,P,a,t)d N qd N p = 1. (6) 

The evolution equation of the distribution function pN(q,p,a,t) is Liouville equation in the 
Hamilton picture 

dp N (q,p,a,t) 

j t = -tt{q,p,a,t)p N {q,p,a,t). (7) 

This equation describes the change of the distribution function p N along the trajectory in the 
6N-dimensional phase space. Here, Q is defined by 

x dFi dGi 

n(q,p,a,t) = -^ + —^. 8 
dpi dqi 

Here and later we mean the sum on the repeated index % from 1 to N. Derivative d/dt is a total 
time derivative 

- = - + G t — + F t —. (9) 

dt dt 1 dqi 1 dpi 

If the vector function Gi is defined by Eq. (4), then 

s dFi d 2 H 
M{q,p,a,t) — — h 



dpi dqidpi ' 

In general, the second term does not vanish, for example, in the nonlinear one-dimensional 
sigma-model [27]. 

In the Liouville picture (for the Lagrange variables) the function Q defines the velocity of 
the phase volume change [26] 



dV ph (a, t) 
dt 



£l(q,p,a,t)d qd p. 



Definition 3. // Q < for all phase space points (q,p) and Q < for some points of phase 
space, then the system is called a dissipative system. 
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We can define dissipative system using a phase density of entropy 

S(q,p,a,t) = —k In pjsr(q,p,a,t). 

This function usually called the Gibbs phase. Eq. (7) leads to the equation for the entropy 
density (Gibbs phase) 

dS(q,p,a,t) 

— = kn(q,p,a,t). (10) 

It is easy to see that the function Q is proportional to the velocity of the phase entropy density 
change. Therefore, the dissipative systems can be defined by the following equivalent definition. 

Definition 4. A system is called a generalized dissipative system if the velocity of the en- 
tropy density change does not equal to zero. 

Let us define the special class of the non-Hamiltonian systems with distribution functions 
that are completely characterized by the Hamiltonian. These distributions can be derived an- 
alytically as stationary solutions of the Liouville equation for the non-Hamiltonian system. 

Definition 5. A non-Hamiltonian system will be called a canonical non-Hamiltonian system 
if the distribution function is determined by the Hamiltonian, i.e., pN(q,p,a) can be written in 
the form 

p N (q,p,a) = p N (H(q,p,a),a), (11) 
where a is a set of external parameters. 

Examples of the canonical non-Hamiltonian systems: 

(1) The constant temperature systems [2, 3, 4, 5, 6, 9] that have the canonical distribution. In 
general, these systems can be defined by the non-holonomic constraint, which is suggested in 
[11]. 

(2) The Fermi-Bose canonical-dissipative systems [13] which are defined by the distribution 
functions in the form 

P»W™> a » = exp{3(H(q,p)a)-p) + s- ^ 

(3) The classical system with the Breit-Wigner distribution function is defined by 

PN(H) = (H-E)* + (172)2 • (13) 

3 Distribution as a Solution of the Liouville Equation 

3.1 Formulation of the Results 

Let us formulate the proposition that allows us to answer the following question: Is this system 
a canonical non-Hamiltonian system? 

Let us consider the N-particle non-Hamiltonian systems which are defined by the equations 

dri dH dpi dH i ^ (n) 



dt dp,' dt dri 



+ F\ n> . (14) 
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The power of non-potential forces is defined by 



,(n) dH 



V(r, P ,a) = Fr — . (15) 

If the power of the non-potential forces is equal to zero (V = 0) and dH/dt = 0, then classical 
system is called a conservative system. The velocity of an elementary phase volume change Q 
is defined by the equation 

n{rpa) = ^l + d " H = ^ (16) 
dpi dridpi dpi 

We use the following notations for the scalar product 

dAj = ",8A^ | dAyj | dA zi 
dsn ~^da xi da yi da zi 

The aim of this section is to prove the following result. 

Proposition 1. If the non-potential forces of the non-Hamiltonian system (14) satisfy 
the constraint condition 

9{H)F * 9p7~^p7-°' (17) 

then this system is a canonical non-Hamiltonian system with the distribution function 

p N (r, p, a) = Z(a)exp[-L(H(r, p, a))], (18) 
where the function L(H) is defined by the equation 



The condition (17) can be formulated in other words: 
// velocity of the elementary phase volume change Q is directly proportional to the power V of 
non-potential forces of the non-Hamiltonian system (14) and coefficient of this proportion- 
ality is a function g(H) of Hamiltonian H, i.e., 

Q(r,p,t)-g(H)V(r,p,t) = 0, (20) 

then this system is a canonical non-Hamiltonian system. 

Note that any non-Hamiltonian system with the non-holonomic constraint (20) or (17) is a 
canonical non-Hamiltonian system. 

Example. Let us consider g(H) = 3N(3(a), where /3(a) = l/kT(a). This case is considered 
in [11]. If we consider the N-particle system with the Hamiltonian 

N D 2 

H(r,p,a)=j:^ + U(r,a), (21) 
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and a linear friction, which is defined by the non-potential forces 



F, (n) = -TP;, (22) 
then the non-holonomic constraint (17) has the form 

E^=*T(«)> (23) 
i=i 

i.e., the kinetic energy of the system must be a constant. The constraint (23) is a non-holonomic 
minimal Gaussian constraint [11, 9]. 

If the function g(H) is defined by g{H) = 3N/3(a), then the non-Hamiltonian system can 
have the canonical Gibbs distribution [11]. The classical systems that are defined by Eqs. (21) 
- (23) are canonical non-Hamiltonian systems. 



3.2 Proof of the Result 

Solving the Liouville equation with the non-holonomic constraint (17), we can obtain the 
(canonical and non-canonical) distributions that are defined by the Hamiltonian. 

Let us consider the Liouville equation for the N-particle distribution function p N = pAr(r, p, a, t). 
This distribution function p N express a probability that a phase space point (r, p) will appear. 
The Liouville equation for this non-Hamiltonian system 

dp N d ( \ d 



expresses the conservation of probability in the phase space. Here, we use 

_ OH OH , n) 

dpi UTi 

We define a total time derivative along the phase space trajectory by 

dt dt ^ dr i % <9p, ^ ^ 

Therefore Eq. (24) can be written in the form (7) 

dpN 



dt 



= -tops, (26) 



where the omega function is defined by Eq. (16). In classical mechanics of Hamiltonian systems 
the right-hand side of the Liouville equation (26) is zero, and the distribution function does not 
change in time. For the non-Hamiltonian systems (14), the omega function (16) does not vanish. 
For this system, the omega function is defined by Eq. (16). For the canonical non-Hamiltonian 
systems, this function is defined by the constraint (17) in the form 



n = g(H)F\ 



(n) dH 



dpi 
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In this case, the Liouville equation has the form 



~df = ~ 9{H)F * w, pN - (27) 



Let us consider the total time derivative of the Hamiltonian. Using equations of motion (14), 
we have 

dH dH dHdH , dH „ M \dH dH „ M dH 

= 1 h ( h F' ) = h F' ' . 

dt dt dpi dvi ^ dri % ' dpi dt 1 dpi 

If dH /dt = 0, then the power V of non-potential forces is equal to the total time derivative of 
the Hamiltonian 

F(n) dH = dH 
% dpi dt 

Eq. (27) can be written in the form 

^-(<- <»> 

Let us consider the following form of this equation: 

If g(H) is an integrable function, then this function can be represented as a derivative 

In this case, we can write Eq. (29) in the form 

din pn dL(H) 



, (31) 
dt dt v ' 

As a result, we have the following solution of the Liouville equation: 

PAr(r, p, a) = Z(a)exp - L(H(t, p, a)). (32) 

The function Z(a) is defined by the normalization condition. It is easy to see that the distri- 
bution function of the non-Hamiltonian system is determined by the Hamiltonian. Therefore, 
this system is a canonical non-Hamiltonian system. 

Note that iV is an arbitrary natural number since we do not use the condition N » 1 or 
TV -> oo. 



4 Non-Holonomic Constraint for Non-Hamiltonian Sys- 
tems 

4.1 Formulation of the Result 

Let us formulate the proposition which allows us to derive the canonical non-Hamiltonian 
systems from any equations of motion of non-Hamiltonian systems. 



8 



The aim of this section is to prove the following result. 



Proposition 2. For any non-Hamiltonian system which is defined by the equation 

dvj _ dH_ dpi_ _ 
dt dp/ dt 

where Fj is the sum of potential and non-potential forces 



dt dpi ' dt 

where the non-potential forces ~F™ ew are defined by 



The vectors Aj and Bj are defined by the equations 



and 



(33) 



F^-^ + F^, (34) 

there exists a canonical non-Hamiltonian system that is defined by the equations 

dri dH dpi 



F™ w , (35) 



F n ew = AfcAfcV- AjAj F _ AjBj dH 
1 A k A k 3 A k A k dpj' 



= dg(H) dH dH F(B) + 8F^dH + 
dH dpi dpj 3 dpi dpj 

+9 (7/)Ff/|- - f^fi, (37) 
dpidpj dpidpj 

<9iJ <9rj 9pj 3 dr,i dpj 

+#)Fr»^-^. (38) 



Note that the forces that are defined by Eqs. (36) - (38) satisfy the non-holonomic constraint 
(20), i.e., 

n(ff)F ne W dH dF T - PH 

9{H)F * dp-, -JJp- J^dV, °- (39) 

4.2 Proof. Part I 

In this subsection, we prove Eq. (36). 

Let us consider the N-particle classical system in the Hamilton picture. Denote the position 
of the ith particle by and its momentum by pi. 

Suppose that the system is subjected to a non-holonomic (non-integrable) constraint in the 
form 

/(r,p) = 0. (40) 
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Differentiation of Eq. (40) with respect to time gives a relation 

A,(r,p)^ + B l (r,p)^ = 0, (41) 

where 

A,<r,p) = J£, B,(r,p) = |. (42) 
An unconstrained motion of the ith particle, where % = 1, N, is described by the equations 

where Fj is a resulting force, which acts on the ith particle. 

The unconstrained motion gives a trajectory which leaves the constraint hypersurface (40). 
The constraint forces Rj must be added to the equation of motion to prevent the deviation 
from the constraint hypersurface 

* =G> , f = F, + R4 . (44) 
The constraint force Rj for the non-holonomic constraint is proportional to the A, [16] 

R = AA l5 (45) 

where the coefficient A of the constraint force term is an undetermined Lagrangian multiplier. 
For the non-holonomic constraint (40), equations of motion (43) are modified as 

= G i} ^ = F, + AAj. (46) 
dt dt K J 

The Lagrangian coefficient A is determined by Eq. (41). Substituting Eq. (44) into Eq. 

(41), we get 

A,(F, + AA,) + BjGj = 0. (47) 
Therefore, the Lagrange multiplier A is equal to 

A,-F,- + B,G, 

A = ~ — * . 48 

A fe A fc 

As a result, we obtain the following equations: 

*2- G - dPi ~ F A AjFj + BjGj (49) 
dt l ' dt ~ 1 1 A k A k ■ [ } 

These equations we can rewrite in the form (43) 

d r i _ n dpi ™ ew ,_„s 



with the new forces 



.„«„. _ AfcAfedy AjAj A , 1 1 j 



F new = — A » J , J F, - -p-^-Gj- (51) 

In general, the forces ~F™ ew are non-potentials forces (see examples in [11]). 

Eqs. (49) are equations of the holonomic non-Hamiltonian system. For any trajectory of 
the system in the phase space, we have / = const. If initial values r fe (0) and Pfc(0) satisfy the 
constraint condition /(r(0),p(0)) = 0, then solution of Eqs. (49) and (51) is a motion of the 
non-holonomic system. 
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4.3 Proof. Part II. 

In this subsection we prove Eqs. (37) and (38). 

Let us consider the non-Hamiltonian system (43) with 

*=-f < 52 > 

and the special form of the non-holonomic constraint (40). Let us assume the following con- 
straint: the velocity of the elementary phase volume change Q(r, p, a) is directly proportional 
to the power V(r, p, a) of the non-potential forces, i.e., 

n(T,p,a)=g(H)V(T,p,a), (53) 

where g(H) depends on the Hamiltonian H. Therefore, the system is subjected to a non- 
holonomic (non-integrable) constraint (40) in the form 

/(r, p, a) = g(H)V(r, p, a) - Q(r, p, a) = 0. (54) 

This constraint is a generalization of the condition which is suggested in [11]. The power V of 
the non-potential forces is defined by Eq. (15). The function Q is defined by Eq. (8). 
Eq. (54) for the non-potential forces has the form 



OH d¥f ] 
d Pj d Pj 



Let us find the functions A, and Bj for this constraint. Differentiation of the function 
f(r, p, a) with respect to p^ gives 



df _ d ( , H)F (n)dH, d gig 

dpi dpi^ 3 dp/ dpi dpj 



Therefore, we obviously have (37). Differentiation of the function f(r,p,t) with respect to r.j 
gives 



B - 9f - 9 (a(H)F^ dH ) 9 ^ 

Bi -wrdFS 9{H)Fj dp-J-fr^ 

Therefore, we have (38). 
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4.4 Minimal Constraint Models 

To realize simulation of the classical systems with canonical and non-canonical distributions, 
we must have the simple constraints. Let us consider the minimal constraint models which are 
defined by the simplest form of the Hamiltonian and the non-potential forces 

if (r, p, a ) = ^- + U(r, a), F\ n} = - 7Pi . (55) 
2m 
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where p 2 = YliLi Pi- For these models, the non-holonomic constraint is defined by the equation 

f = g(H)^-3N = 0, (56) 

where N is the number of particles. The phase space gradients (37) and (38) of the constraint 
can be represented in the form 

The non-potential forces of the minimal constraint models have the form 

w _ P 2 $ij ~ ViVj dU PiPj dg(H) dU 



p 2 drj 2p 2 {(p 2 /2m){dg{H)/dH) + g(H)) dH drj' 

It is easy to see that all minimal constraint models have the potential forces. 
For the minimal Gaussian constraint model 

dg(H) 

OH 

we have the non-potential forces in the form 

dU p 2 5ij - PiPj 



F; = 



drj p 2 



This model describes the constant temperature systems [2, 3, 4, 5, 6, 9, 11]. 

4.5 Minimal Gaussian Constraint Model 

Let us consider the N-particle system with the Hamiltonian 

H(r,p,a) = ^- + U(r,a), (57) 
the function g(H) = 3N/kT, and the linear friction 

F\ n) = -TP;, (58) 

where % = 1,...,N. Note that is an arbitrary natural number. Substituting Eq. (58) into 
Eqs. (15) and (16), we get the power V and the omega function fi: 

V = -— p 2 , n = -3-fN. 
m 



The non-holonomic constraint has the form 

^2 



— = kT(a), (59) 
m 
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i.e., the kinetic energy of the system must be a constant. Note that Eq. (59) has not the 
friction parameter 7. 

For the N-particle system with friction (58) and non-holonomic constraint (59), we have the 
following equations of motion 

dri pi dpi dU df 
dt m dt OYi dpi 

where the function / is defined by 

f(v,p)= l -{p 2 -mkT) : /(r,p) = 0. (61) 

Eq. (60) and condition (61) define 6N+1 variables (r, p, A). 

Let us find the Lagrange multiplier A. Substituting Eq. (61) into Eq. (60), we get 

dpi dU s , nnS 



Using df /dt = in the form 



Pi f = (63) 



and substituting Eq. (62) into Eq. (63), we get the Lagrange multiplier A in the form 

1 dU 

As a result, we have the holonomic system that is defined by the equations 



cfci _ pi dpi _ 1 dU_ _dU_ . , 

dt ~ m dt ~ mkT PiPj drj dr t ' 1 j 

This system is equivalent to the non-holonomic system (60). For the classical N-particle system 
(64), condition (59) is satisfied. If the time evolution of the N-particle system is defined by Eq. 
(64) or Eqs. (60) and (61), then we have the canonical distribution function in the form 

p(r, p, a, T) = exp^ (>(a, T) - if (r, p, a)) . (65) 

For example, the N-particle system with the forces 

uj 2 (a) 9 . 
F; = —j-^-piPjTj - mu (a)ri (66) 

can have canonical distribution (65) of the linear harmonic oscillator with 



U(r,a) 



mu 2 (a)v 2 
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5 Canonical Distributions 



In this section, we consider the subclass of the canonical non-Hamiltonian system that is de- 
scribed by canonical distribution. This subclass of the canonical non-Hamiltonian N-particle 
system is defined by the simple function g(H) = 3Nf3(a) in the non-holonomic constraint (20). 

Proposition 3. If velocity of the elementary phase volume change is directly proportional 
to the power of non-potential forces, then we have the usual canonical Gibbs distribution as a 
solution of the Liouville equation. 

In other words, the non-Hamiltonian system with the non-holonomic constraint 

n = (3(a)V (67) 

can have the canonical Gibbs distribution 

Pn = exp(3(a) [J 7 (a) — H(r, p, a) 



as a solution of the Liouville equation. Here, the coefficient (3(a) does not depend on (r, p, t), 
i.e., 

d(3(a)/dt = 0. 

For the non-Hamiltonian systems, the omega function (16) does not vanish. Using Eq. (15), 
we have 

n = /?(a)F< B >^. (68) 
opt 

In this case, the Liouville equation has the form 



Let us consider the total time derivative for the Hamiltonian 

dH dH „ M dH 

ir = flr +B *V (70) 

If dH /dt = 0, then the energy change is equal to the power V of the non-potential forces F^. 
Eq. (69) can be written in the form 

dpN n, \dH . 
If = -M-dt PN - (71) 

Therefore, the Liouville equation can be rewritten in the form 

d lnp N (T,p,a,t) dH(r, p, a) = 

dt [a) dt 

Since coefficient (3(a) is a constant (d(3(a)/dt = 0), we have 

j t (lnp N (r, p, a, t) + (3(a)H(r, p, a)) = 0, 

14 



i.e., the value (InpN + (3H) is a constant along the trajectory of the system in 6N-dimensional 
phase space. Let us denote this constant value by /3(a) J 7 (a). Then, we have 

lnp N (r, p, a, t) + @(a)H(r, p, a) = f3(a)J r (a), 
where d!F(a)/dt = 0. It follows that 

In pAr(r, p, a, t) = (3(a) (^F(a) — H(r, p, a)) . 
As a result, we have a canonical distribution function 

p N (r, p, a, t) = exp(3(a) (.F(a) - H(r, p, a)) 

in the Hamilton picture. The value T(a) is defined by the normalization condition (40). 

Therefore the distribution of this non-Hamiltonian system is a canonical distribution. 

Note that N is an arbitrary natural number since we do not use the condition N » 1 or 
TV -> oo. 



6 Non-Canonical Distributions 

The well-known non-Gaussian distribution is the Breit-Wigner distribution. This distribution 
has a probability density function in the form 

The Breit-Wigner distribution is also known in statistics as Cauchy distribution. The Breit- 
Wigner distribution is a generalized form originally introduced [28] to describe the cross-section 
of resonant nuclear scattering in the form 

P{H) = (H-EY + {T/2f (73) 

This distribution can be derived from the transition probability of a resonant state with known 
lifetime [29, 30, 31]. 

If the function g(H) of the non-holonomic constrain is defined by 

*w - ( H 2 -"i + E iTi2r (74) 

then we have non-Hamiltonian systems with the Breit-Wigner distribution as a solution of the 
Liouville equation. 

If the function g(H) of the non-holonomic constrain has the form 

9(H) = l + aexp(3(a)H> (75) 

then we have classical non-Hamiltonian systems with Fermi-Bose distribution (12) considered 
by Ebeling [13]. This distribution can be derived as a solution of the Liouville equation. Note 
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that Ebeling derives the Fermi-Bose distribution function as a solution of the Fokker-Planck 
equation. It is known that Fokker-Planck equation can be derived from the Liouville equation 
[14]- 

If the non-potential forces F\ n ^ are determined by the Hamiltonian 

F< n) = -dG(H)/d Pt , (76) 

then we have the canonical non-Hamiltonian systems, which are considered in [12, 13]. These 
systems are called canonical dissipative systems. 
Note that the linear function g(H) in the form 

g{H) = A (a) + fa(a)H 

leads to the following non-canonical distribution function 

p N = Z(a)exp - (A(a)if + ^f3 2 (a)H 2 ). (77) 

The proof of this proposition can be directly derived from Eqs. (32) and (30). 
Let us assume that Eq. (11) can be solved in the form 

H = 0(a)h(p N ), (78) 

where h depends on the distribution p N . The function 6(a) is a function of the parameters a. 
In this case, the function g(H) is a composite function 

R(p N ) = -g(9(a)h(p N )). (79) 

This function can be defined by 

R(pN) = ±( H a) d -!^.y'. (80) 
Pn \ op N J 

In this case, the Liouville equation for the non-Hamiltonian system has the form 

dpN 



dt 



R{p N )V. (81) 



This equation is a nonlinear equation. Note that the classical Fermi-Bose systems [13] have the 
function in the form 

R(p N ) = -P(a)(p N -sp 2 N ). (82) 

The nonlinearity of the Liouville equation is not connected with an incorrectly defined 
phase space. This nonlinearity is a symptom of the use of an incorrectly defined boundary 
condition. The Bogoliubov principle of correlation weakening cannot be used for classical 
Fermi-Bose systems. The classical Fermi-Bose systems can be considered as a model of open 
(non-Hamiltonian) system with the special correlation. Note that the nonlinear evolution of 
statistical systems is considered in [33, 34, 35, 36, 37, 38, 39]. 
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7 Thermodynamics Laws for non-Hamiltonian Systems 

Let us define the mean value f(a) of the function /(r, p, a) by the relation 

f( a ) = J /( r > P> a)pN(r, p, a)rf Jv rrf 7V p (83) 
and the variation for this function by 

S a f(r,p,a) = j2 df{ l^ a) da k . (84) 

The first law of thermodynamic states that the internal energy U(a) may change because 
of heat transfer SQ, and work of thermodynamics forces 

n 

5A = Y J F k (a)da k . (85) 
k=i 

The external parameters a = {ai,a 2 , ...a n } here act as generalized coordinates. In the usual 
equilibrium thermodynamics the work done does not entirely account for the change in the 
internal energy. The internal energy also changes because of the transfer of heat, and so 

dU = 5Q- 5 A. (86) 

Since thermodynamics forces F k (a) are non-potential forces 

dF k (a) 9F,(a) 



dai da k 



(87) 



the amount of work 5 A depends on the path of transition from one state in parameters space 
to another. For this reason 5 A and SQ, taken separately, are not total differentials. 

Let us give a statistical definition of thermodynamic forces for the non-Hamiltonian systems 
in the mathematical expression of the analog of the first thermodynamics law for the mean 
values. It would be natural to define the internal energy as the mean value of Hamiltonian 

17(a) = J H(r,p,a)p N (r,p,a)d N rd N p. (88) 

It follows that the expression for the total differential has the form 

dU(a) = J S a H(r,p,a)p N (r,p,a)d N rd N p + J H(r,p,a)S a p N (r,p,a)d N rd N p. 

Therefore 

dU(a) = J dH{r d ^ a) Sa kPN {v, p, a)d N vd N p + J H(r, p, a)S a p N {v, p, a)d N vd N p. (89) 

In the first term on the right-hand side, we can use the definition of phase density of the 
thermodynamics force 

?phu n ^ ^ dH(r,p,a) 
da k 
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k (r,P,a) 



The thermodynamics force Fk(a) is a mean value of the phase density of the thermodynamics 
force 

F k (a) = J F p k h (r, p, a)p N (v, p, a)d N vd N p. (90) 

Using this equation we can prove relation (87). 

Analyzing these expressions we see that the first term on the right-hand side of Eq. (89) 
answers for the work (85) of thermodynamics forces (90), whereas the amount of the heat 
transfer is given by 

5Q = J H(r, p, a)5 aPN (r, p, a)d N vd N p. (91) 

We see that the heat transfer term accounts for the change in the internal energy due not to 
the work of thermodynamics forces, but rather to change in the distribution function cased by 
the external parameters a. 

Now let us turn our attention to the analog of the second law for the non-Hamiltonian 
systems. 

The second law of thermodynamics has the form 

5Q = 6(a)dS(a). (92) 

This implies that there exists a function of state S(a) called entropy. The function 9(a) acts as 
integration factor. 

Let us prove that (92) follows from the statistical definition of 5Q in Eq. (91). For Eq. 
(91), we take the distribution that is defined by the Hamiltonian, and show that (91) can be 
reduced to (92). 

Let us assume that Eq. (11) can be solved in the form 

H = 9(a)h(p N ), 

where h depends on the distribution p^. The function 9(a) is a function of the parameters 
a = {ai, 02, a n }. 

We rewrite (91) in the equivalent form 

5Q = J (9(a)h(p(r, p, a), a) + C(a))5 a p N (r, p, a)d N rd N p. (93) 

New term with C (a) , which is added into this equation, is equal to zero because of the normal- 
ization condition of the distribution function p^: 

C(a)5 a j p N (r,p,a)d N rd N p = C(a)5 a l = 0. 

We can write Eq. (93) in the form 

5Q = 9(a)5 a J K(p N (r, p, a))d N rd N p, (94) 

where the function K = K(p^) is defined by 

^^l = h(p N ) + C(a)/9(a). (95) 
op N 



We see that the expression for 5Q is integrable. If we take l/9(a) for the integration factor, 
thus identifying 9(a) with the analog of absolute temperature, then, using (92) and (94), we 
can give the statistical definition of entropy 

S(a) = J K{p N {r, p, a))d N rd N p + S . (96) 

Here, So is the contribution to the entropy which does not depend on the variables a, but may 
depend on the number of particles iV in the system. Note that the expression for entropy is 
equivalent to the mean value of phase density function 

S ph (r, p, a) = K(p N (r, p, a))/p N {r, p, a) + C(a). (97) 

S ph is a function of dynamic variables r, p, and the parameters a = {a 1? a2, a n }. The number 
N is an arbitrary natural number since we do not use the condition N >> 1 or N — > oo. Note 
that in the usual equilibrium thermodynamics the function 9(a) is a mean value of kinetic 
energy. In the suggested thermodynamics for the non-Hamiltonian systems 9(a) is the usual 
function of the external parameters a = {ai, 02, a n }. 

8 Conclusion 

The aim of this paper is the extension of the statistical mechanics of conservative Hamiltonian 
systems to non-Hamiltonian and dissipative systems. In this paper, we consider a wide class of 
non-Hamiltonian statistical systems that have (canonical or non-canonical) distributions that 
are defined by Hamiltonian. This class can be described by the non-holonomic (non-integrable) 
constraint: the velocity of the elementary phase volume change is directly proportional to the 
power of non-potential forces. The coefficient of this proportionality is defined by Hamiltonian. 
The special constraint allows us to derive solution for the distribution function of the system, 
even in far-from equilibrium situation. These distributions, which are defined by Hamiltonian, 
can be derived analytically as solutions of the Liouville equation for non-Hamiltonian systems. 

The suggested class of the non-Hamiltonian systems is characterized by the distribution 
functions that are determined by the Hamiltonian. The constant temperature systems [2, 3, 
4, 5, 6, 9], the canonical-dissipative systems [12, 13], and the Fermi-Bose classical systems [13] 
are the special cases of suggested class of non-Hamiltonian systems. For the non-Hamiltonian 
N-particle systems of this class, we can use the analogs of the usual thermodynamics laws. Note 
that iV is an arbitrary natural number since we do not use the condition N » 1 or N — > 00. 
This allows us to use the suggested class of non-Hamiltonian systems for the simulation schemes 
[32] for the molecular dynamics. 

In the papers [40, 41, 42], the quantization of the evolution equations for non-Hamiltonian 
and dissipative systems was suggested. Using this quantization it is easy to derive the quantum 
Liouville- von Neumann equations for the N-particle statistical operator of the non-Hamiltonian 
quantum system [26]. We can derive the canonical and non-canonical statistical operators that 
are determined by the Hamiltonian [17, 18]. The condition for non-Hamiltonian systems can 
be generalized by the quantization method suggested in [40, 41]. 
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